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Abstract 



Vh ' Covariant equations characterizing the strength of a singularity in spherical 

symmetry are derived and several models are investigated. The difference 

between central and non-central singularities is emphasised. A slight modifi- 

f^ , cation to the definition of singularity strength is suggested. The gravitational 

^vq \ weakness of shell crossing singularities in collapsing spherical dust is proven 

K^ ' for timelike geodesies, closing a gap in the proof. 
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I. INTRODUCTION 

Over thirty years have passed since the appearance of the first of a series of theorems 
estabhshing that under very general circumstances, space-times must develop singularities. 
This first result, due to Penrose |jl|, appeared in 1965 and the body of work which grew 
up around the singularity theorems is contained in the book of Hawking and Ellis [§ which 
was first published in 1973. However, out understanding of the nature of these singularities 
remains far from complete. This lack is best exemplified by the absence of a proof, or 
definitive refutation, of the cosmic censorship conjecture (CCC) 0,^. 

An important aspect of a singularity is its gravitational strength ||^. A singularity is 
termed gravitationally strong, or simply strong, if it destroys by crushing or stretching any 
object which falls into it. The most familiar example is the singularity at r = in the 
Schwarzschild solution. (Throughout this paper, we will refer to a singularity at r = as 
a central singularity, and to others as non-central.) A radially inf ailing object is infinitely 
stretched in the radial direction and crushed in the tangential directions, with the net result 
of crushing to zero volume. A singularity is termed weak if no object which falls into the 
singularity is destroyed in this way. The mathematical description of these ideas runs as 
follows [|,|. 

Let 7 : [to,0) — i> M be a causal geodesic which approaches a singularity as r — > 0^. 
Define Juil) for u G [to,0) to be the set of maps Z : [z^, 0) -^ TM satisfying the geodesic 
deviation equation along 7 such that Z{t) G T^{t-)M, gabk"'{T)Z^{T) = where fc" is the 
tangent to 7 and 



M = o 



i.e. Juil) is the set of Jacobi fields along 7 which vanish at 7(2^). Along a timelike geodesic, 
three independent Jacobi fields define, via the exterior product, a volume element V{t) along 
7. Along a null geodesic, two such fields define an area element which we also denote V{t). 
A timelike (null) geodesic terminates in a strong curvature singularity if for all u G [tq, 0) 
and all independent triads (dyads) in Jui'j) we have 

lim_ V{t) =0. 

Then the singularity itself is said to be strong if every causal geodesic which approaches it 
terminates in a strong singularity. The geodesic terminates in a weak singularity if the limit 
above is finite and non-zero, and the singularity is weak if every causal geodesic approaching 
it terminates weakly. We will argue below for a slight modification of this definition whereby 
the term strong is also attached to a singularity if the norms of the Jacobi fields themselves 
have zero or infinite limit. 

The importance of the notion of the gravitational strength of a singularity for the CCC 
is that a statement of such possibly need not rule out the occurrence of naked weak singular- 
ities. This is based on the belief that one may extend the geometry of space-time through a 
weak singularity without traumatic effects f^J^. A general description of this extension does 
not exist - indeed as far as the author can determine, only two examples of this procedure 
exist in the literature, one due to Papapetrou and Hamoui |Q and the other due to Clarke 
and O'Donnell |p. Both deal with extending through a shell crossing singularity in collaps- 
ing spherical dust. However the fact that at a weak singularity one has, along any timelike 



geodesic, a finite non-degenerate triad of Jacobi fields, from wliicli it may be possible to 
construct a metric in a canonical way, lends support to the idea. 

Our aim here is to give a comprehensive analysis of the strengths of singularities in 
spherical symmetry. This has been the arena of some of the most interesting developments 
in general relativity in recent years, and an understanding of what can and cannot occur 
in spherical symmetry may be a valuable guide for more general situations. Specifically, 
we study the Jacobi equations for arbitrary radial causal geodesies. This allows us to give 
covariant equations identifying the geometrical terms which control the strength of the 
singularity. As one would expect, the results are simpler than those obtained by Clarke 
and Krolak |T0| which apply to general space-times. We study three different models which 



help illustrate the different situations which occur, and by way of application, demonstrate 
that (i) a non-central singularity is always weak along null directions and (ii) the shell- 
crossing singularities in collapsing spherical dust are weak. (This has only been demonstrated 
previously for radial null directions |TT|; we complete the proof by showing that it is also 
true for radial timelike directions.) 



II. RADIAL CAUSAL GEODESICS, JACOBI FIELDS AND THE VOLUME 

ELEMENT. 

The line element of a spherically symmetric space-time may be written as 

ds^ = -2e~^^du dv + r'^du'^, (2.1) 

where / = /(u, f ), r = r{u, v) and duj"^ is the line element of the unit 2-sphere. The function 
r is an invariant of the space-time which we will, quite properly, refer to as the radius. 
The coordinates u and v are both null, labelling the null hypersurfaces generated by the 
two families of null geodesies orthogonal to the orbits of the 5*0(3) symmetry group of the 
space-time. The form ( p.l|) is invariant under the transformations u — > u'{u), v -^ v'{v). 
The non-vanishing Ricci tensor components are 

(2.2a) 
(2.2b) 
(2.2c) 
(2.2d) 

We use the convention that subscripts attached to lower case letters refer to partial 
derivatives, but elsewhere refer to tensor components in the associated coordinate basis. 
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The Misner-Sharp energy |T2| is 



E = ^(l + 2e2/r„r,), (2.3) 

and the Weyl tensor is completely determined by the Newman- Penrose term \l/2, calculated 
on a null tetrad based on the principal null directions of the space time. Thus \l/2 is an 
invariant of the space-time and is given by 

^2 = |^(r™ + rA„)-^. (2.4) 



We note the further invariant, 



Of „ -C/ -ft , . 

e'^/™ = - + 2^2--j, 2.5 

r'^ 12 

where R is the Ricci scalar. 

The radial geodesic equations are 

- 2e-^^uv = e (2.6a) 

u - 2fuV? = (2.6b) 

V - 2fy = (2.6c) 

where e = for null geodesies and —1 for timelike; the overdot is respectively, differentiation 
with respect to afiine parameter and proper time. Space-like geodesies will not concern us 
here. 

We now look at the Jacobi fields along arbitrary radial causal geodesies, beginning with 
time-like geodesies. 

The unit tangent to an arbitrary time-like curve 7(r) in the radial 2-space can be written 
in the form 

- , d 1 2f, -1 ^ 

^ = ^TT + o^^ V' 
ou 2 ov 

where h = h{u,v). The condition that 7 be geodesic is then 

K + 2h\-^f{K-2Uh) = Q. (2.7) 



This follows from the geodesic equations ( ^.61) . The variation of any scalar quantity s along 
this geodesic is given by 



and using 



s = k^VaS = hsu + -e^^h-h^, (2.8) 



S — n ySuu ~r '^JuSuJ ' A^ '^ ySvv ~r ^Jv^vj ~r C S^v \^''^) 



A Jacobi field Z"" along 7 satisfies the geodesic deviation equation 

Z'' + RM''Z^k^k'^ = Q, (2.10) 

which is a linear equation for Z"' and so a basis for the Jacobi fields may be found by 
obtaining all independent Jacobi fields in the radial 2-space and in the tangential 2-space. 
We take 

- (9 - d 

^(1) =x(M,t;)— ii2) = y{u,v) CSC e— 

as candidates for the Jacobi fields in the tangential 2-space. Note that the norms of ^[\) and 
^?2) are |x| and \y\ respectively. The geodesic deviation equation ( p.lO|) applied to ^?^x yields 
the following equation for x (the same result applies to y): 






Using (2^) and (|2.9| ), this assumes the remarkably simple form 

rx + 2rx = 0, (2.11) 

which can be integrated to obtain 

where Xq is constant and we have included the initial condition x(ri) = 0, so that ^(i)(ti) = 0. 
The second linearly independent solution is x(r) = 0. We obtain the same result for ,^(2): 

Before dealing with the radial Jacobi fields along the time-like geodesies, we describe the 
situation for radial null geodesies. It turns out to be remarkably simple. From (p. 6a]) , either 



u 01 V vanishes along a radial null geodesic. Take it to be the latter. We can then integrate 
( p.6b|) to obtain ii = ce^^ and so the tangent is 



r 2f ^ 

k = ce ■' —-. 

ou 
The variation of a scalar s along the geodesic is s = ce^^ Su and 

In the null case, there are only tangential Jacobi fields. Candidates for such are given by 
^(i_2) as above, and it turns out that the norms x, y obey the same equation as in the time-like 
case, and thus the solutions are given by ( |2.12a|) . 

We now treat the radial Jacobi fields along a time-like geodesic. A space-like vector in 
the radial 2-space orthogonal to k"- has the form 

-> . d I Of J ^1 d 

t, = an— e ■' ah -—, 

ou 2 ov 

where a = a{u, v). This has norm \a\. Using ( p.9|) , we find that the condition for ^ to satisfy 
the geodesic deviation equation along k"- is 

a + 2e2^/™a = 0. (2.13) 

Recall that according to (^^.5]), e^^ fuv is an invariant of the space-time. Thus equations 
( p.l2| ) and ( p.l3|) provide a covariant set of equations which will determine the strength of 
the singularity. To see how this comes about, we obtain the relationship between V{t) and 
the quantities a, x, y. 



Since each of a, x and y satisfy second order linear ordinary differential equations, there 
are six independent Jacobi fields along a time-like geodesic 7. An arbitrary triad of corre- 
sponding l-forms is given by 



Zq = ttaG + Xa.r'^dO + y^ sin 

where e = {2h)^^du — e^'^^hdv and a = 1,2,3. Then the general volume element along 7 
has the form 

V{r) = zi Az2 Az3 

= 6a[iX2y3]r'^ sin 9e Ad9 A dcp. 

The norm ||W|| of a p— form W = iy[j^...j ]dx^^ A ... A dx^^ is given by 

||Wf = Wh^ i \W'^-'^, 

W W I *1 ••■ tp I J 

where the vertical bars indicate summation only over ii < 12 < ... < ip. This gives 

\\V{r)\\=Q\a[iX2y3y. (2.14) 

The existence of six independent solutions of the geodesic deviation equation indicates 
that in a general space-time, the different V{t) form a six parameter family. This surfeit 
of possibilities would produce problems if one wanted to make statements about singularity 
strengths based on the behaviour of all such V{t). However the definition of J^-^ reduces the 
number significantly. Note that Zq, vanishes if and only if each of a^, x^ and y^ vanish. The 
general solution for a has the form 

a = c+aj^ + C-a^ 



where c± are arbitrary constants and a± are any two independent solutions of ( |2.13|) . A 



similar result holds for x and y. The initial condition a{Ti) = fixes the ratio c_|_/c_, so 
that for the problem in hand, there is only one choice, up to a constant multiple, for each 
of a, X and y. Therefore the norm of every relevant volume element has the simple form 

\\V{T)\\ = \axyy (2.15) 



where a, x and y here represent the general solutions of the appropriate equations ( p.l2|) and 
( p.l3| ), into which constants may be absorbed. 

We thus have a simple and direct way of assessing the strength of a singularity. We deter- 
mine the limiting behaviour of solutions of (|2.12|) and ( p.l3|) as the singularity is approached 
and then use ( |2.15| ) to calculate lim^^Q- 11^(^)11- 



III. POSSIBLE SOLUTIONS OF THE MAIN EQUATIONS 

We now consider the various possible limiting behaviours which can occur for solutions 
of the main equations, i.e. ( |2.12| ) and ( ^^.131 ). 

From ( p.l5| ), the forms rx and ry arise naturally and occur in the following proposition 



which applies to both x and y. 
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Proposition One (i) For a non-central singularity, IIiiIt-^o- "^^ is finite and non-zero, 
(a) For a central singularity, if 



then lim^^o- rx = 0. Otherwise 



lim / r ^(r') dr' < cxd, 



lim rx = — lim 



r^O- r^O- r{T) 



Here and below, the limit refers to the limit as r ^ 0^ along a geodesic which approaches 
the singularity at r = 0. The proofs of part (i) and the first part of part (ii) follow im- 
mediately from ( p.l2a| ), and that of the second part from ( 2.12a ) and I'Hopital's rule. A 



consequence of this is that the strength of a non-central singularity is completely determined 
by the limiting behaviour of a at the singularity. 

Next we summarise the possible behaviour of a in the appropriate limit. Define F{t) = 
2e^^fuvir) and G(r) = t^F{t). Then j ^A^ is equivalent to 



T^'d + G{r)a = 0, (3.1) 

and we wish to determine the behaviour of a at r = which is a singular point of ( |3.1| ). We 



quote the following results from Bender and Orszag [T^, which should also be found in any 
text on linear differential equations. All the asymptotic relations below hold as r — > 0~ and 
c, c± are arbitrary constants. 

The equation is said to have a regular singular point at r = if G{t) is analytic in a 
neighbourhood of r = 0. Otherwise, r = is called an irregular singular point. For a regular 
singular point, we define Gq = G{0). In fact the results below apply more generally, namely 
if G{t) = 0(1) as r — >• 0~ with a very low degree of differentiability; G G C^(— ro,0] is 
sufficient. The method of Frobenius applies. The roots of the indicial equation are 

vi,2 = l±l{i-^Goy/'. 

The following possibilities arise. 

(RSPl) f 1 — f 2 ^ Z. Then a^r) ~ c+|r|''^ + c_|r|''^. Three subcases arise depending on the 

value of Go- 

(RSPla) 1/4 < Go. Then a{T) ~ c|r|i/^ so that lim^_o- a^r) = 0.. 

(RSPlb) < Go < 1/4. Then fi 2 are both positive so that limT-_^o- ci{T) = 0. 

(RSPlc) Go < 0. Then t>2 < so that lim^^o- ^^(t) = cxd. 

(RSP2)t;i = v2^Go = 1/4 Then a(r) ~ c+|r|^/2^c_|r|i/2 in |r|. Again, lim^^o- a(r) = 0. 

(RSP3) vi - V2 e N - ^ Go = {1 - k'^)/4:, k e N+. Then a^r) ~ c+|r|^i + c.{\r\''^ + 

d\T\^^ In \t\) where d is a fixed constant. We mention under this last heading one special case 

of particular importance, that for which Go = 0. This includes singularities whereat F{t) 

is finite and, typically, space-times with weak non-central singularities. 

(RSP3a) Go = 0. Then a{T) ~ c_ -(- c+|r| -f c_c?|r| In |r|), and so lim^^o- ct(T) is finite and 

non-zero. 



The second class of possibilities arises when ( p.l|) has an irregular singular point at r = 0. 
If liniT-^o- |G'(t)| = oo, then the WKB approximation holds. This gives 



a{T) ~ c(F(r))-i/^exp{± [ {F{T')f'^ dr']. 



There are two possibilities here. 

(ISPl) liniT-^o- F{t') = +00. Then liniT-^o- o(''") = +00. 

(ISP2) lim^^o- Fir) = -00. Then lim^^o- air) = 0. 

This does not cover all possibilities since there are irregular singular points for which the 
limit limT-^o- |C^(''")I does not exist. Typically, this would occur if G{t) is oscillatory in a 
neighbourhood of r = 0, e.g. 

Gp{T)=k\T\-'smi\T\-P), 

where p > 0. Taking this form for G{r) and defining x = \t\~^ , b = x^^^^ ^I'^a ( |3.1|) becomes 



6" + 



k sin X 1 — p "^ 

_ — I tl — 

p2 J.2-P 1 42;2 



0. 



The dominant coefficient of h is the decaying oscillatory term, and this determines the 
asymptotic behaviour of the solutions. There are three different cases, depending on the 
value of p. We quote the result for the simplest case, which is p > 1. The asymptotic 
behaviour in this case is given by |T4 



ai ~ x(i-P")/2^ a2 ~ x-(i+P")/2 (x ^ 00). 

Thus the singularity is strong and stretching for p > 1. It turns out that the singularity 



is strong and crushing for 1/2 < p < 1, see JT^ for details. On the other hand, if we take 
G{t) = A;sin(|r|^^|), the same procedure leads to 

,,, fc sinx, 

pZ ^z 



the asymptotic solutions of which lead to 1 14 



ai ~ 1, 02 ~ X 



X — > 00 



and so in this case the singularity is weak. Notice that we have in this case G{t) = 0(1), 
but G is not differentiable at r = 0. 

Our main point here is that both strong and weak singularities may occur in this class 
and the analysis to determine which case obtains may be quite difficult. 
(ISPS) lim^_,o- G{t) does not exist. The singularity may be either strong or weak. 

Keep in mind that the behaviour described here is characteristic of a particular radial 
timelike geodesic which runs into the singularity, and not of the singularity itself. We will 
therefore refer to, for example, a type (RSPla) geodesic, and to a type (RSPla) singularity 
only if all the radial timelike geodesies terminating there are type (RSPla). 

In this language, the central singularity of Schwarzschild space-time is type (RSPlc), 
with a{T) ~ c_|_|r|^/^ + C-\t\'^^^^. Also, rx{T) ~ xoItI-"-/^, ryij) ~ |/o|t|"^^^, so that overall, 
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||l^(r)|| ~ (i|r|^/^, giving a singularity which is strong along timelike approaches. Suppose 
instead the behaviour was rx^r) ~ a;o|r|^/^, ry^r) ~ yokl^''^- Then ||V^(t)|| ~ d (constant), 
so by the current definition, the singularity is weak along tinielike approaches. It would be 
of very little comfort to an observer jumping into such a singularity to realise, as he watched 
his legs elongate and disintegrate, that such volume forms were preserved on his journey. 
The possibility of the existence of such a singularity was noted by Tipler 0. We give an 
example of such below. This motivates the following addendum to the definition of a strong 
singularity. 

We will say that a causal geodesic j : [tq,0) —>■ M approaching a singularity as r ^ 0~ 
terminates in a strong singularity if for all ti G [tq, 0), except some suitably small set (finite, 
countable, zero- measure), the general element of Jnil) is degenerate or infinite in the limit 
T —>■ 0~. We will say that 7 terminates in a weak singularity if the general element of Jnil) 
is finite and non-degenerate in the limit. The terms will be applied to the singularity itself 
if all causal geodesies approaching the singularity behave in one of the two ways. 

By degenerate, we mean that both of the independent Jacobi fields in some particular 
direction (or mutually orthogonal directions) orthogonal to k"" shrink to zero magnitude. A 
non-central type (RSPlb) singularity would be an example of such. 

We now gather the results above into some general statements. 

Proposition Two For a non-central singularity and for a central singularity for which 
r has a finite and non-zero limit along every causal geodesic approaching the singularity, the 
strength of the singularity is determined by i\2.13i ). If the singularity is of type (RSPSa), 
then it is weak. The only other type which may he weak is (ISPS); singularities of the other 
types are strong. A central singularity for which r has zero or infinite limit along every 
causal geodesic approaching the singularity is strong. 

The proof follows from Proposition One and from the definitions above; essentially it 
amounts to some useful book-keeping. A great many singularities will have different be- 
haviours along different geodesies approaching the singularity, and so will not be covered 
by this result. There remains the problem of determining the behaviour of a{T) and r{T) in 
the limit as the singularity is approached. However we have identified which elements of the 
geometry determine the strength of a singularity and listed the various possibilities. 

We now give some applications of the results laid out above. 

IV. EXAMPLES 

In this section, we study the strengths of some singularities in four different (classes) 



of space-times. The first three, two toy models and Roberts' space-time ||T5[, are used to 
illustrate the types of singularities which may arise and some of the points made above. 
The fourth is the marginally bound case of a Lemaitre-Tolman-Bondi (LTB) collapsing dust 



sphere [^. We use the theory above to demonstrate conclusively the weakness of shell- 



crossing singularities in this space-time. 

A. A toy model 
We consider the space-time with line-element 



,2 , , /^-«\^" . 2 

as = —dudv + I — - — I au) , 

i.e. 2e~^-^ = 1 and r = {{v — m)/2)". We take a > 1. The case a = 1 is flat space-time. The 
Ricci scalar is 

i?=4(l-(3«'-2a)^'"'/"), 

and so there is a scalar curvature singularity at r = 0. Since / is constant, the radial Jacobi 
field orthogonal to an arbitrary timelike geodesic will satisfy, according to ( |2.13D , a = 0, 
with general solution a = c+ + C-T. Thus the strength of this central singularity will be 
determined by the behaviour of the tangential Jacobi fields. 

Along a radial null geodesic, we have (without loss of generality) v =constant and u = 
CT + d, where c, d are constants and r is an affine parameter. Thus after a reparametrisation 
of r, we have 

r = clrl". 



The same result holds for all radial timelike geodesies. From ( |2.12| ) and ( |2.15| ), we find 



||y(r)||oc(c+ + c_ 



r T 



\2~2a 



Thus the singularity is strong. Notice that ||^(t)|| ^ cxd as the singularity is approached 
along any radial causal geodesic. The deformation results from infinite tangential stretching. 
The purpose of examining this model is to gave an explicit example where the behaviour 
at the singularity is clearly pathological and destructive, but which would not previously 
have been described as a strong singularity. We note that 

so that the weak energy condition is violated for the values of a of interest here. 



B. Roberts' solution 



Roberts' solution has been used in studies of cosmic censorship |T^ and critical collapse 
||17|| . The line-element is 

ds^ = —dudv + -{u^ — 2uv + (1 — p^)v^)duP'^ 
where p is constant, p = gives fiat space-time. The Ricci scalar is 

^ = ^P ^^' 

and so there is a central scalar curvature singularity. As above, the strength of the singularity 
is determined by the tangential Jacobi fields. In this case we find that r = c\t\ along any 
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radial causal geodesic terminating at r = at parameter value r = 0. We use ( p.l2| ) and 
( p.l5| ) to obtain 



||\/(r)||~A;(c+ + c_|r|) 

with rx, ry ^constant as r ^ 0^, and so this central singularity is weak. Thus the examples 
where this singularity is naked may not be genuine counterexamples to cosmic censorship. 
See also [HI for related comments. 



C. Another toy model 

A model with slightly more complicated dynamics and which illustrates well some of the 
points made above is that with the line element 

as = — \ dudv + duj . 

\ 2 ) \ 2 ) 

We take a > 0; a = is flat space-time. We find that 

2^ffuv = -4«(2a + l)r-^°-^ (4.1) 

and so there is a scalar curvature singularity at r = (recall that this term is an invariant). 
For an arbitrary radial null geodesic, we make take u = Mo=constant. Then we find 

v-u = v-uo = (cr + rf)i/(2a+i)^ 
so that 

after an appropriate shift in the afiine parameter r. Applying the second part of Proposition 
One, we see that all radial null geodesies approaching r = terminate in a strong curvature 
singularity with the area element obeying ||V(r)|| — ;> in every case. 

To solve for the radial timelike geodesies, we make the change of variables r = (f — u)/2, 
t = {v + u)/2. Then the line element takes the form 

ds" = r^^'^-dt' + dr^) + r''du\ 

The geodesic equations for radial infall yield 

r = -r-2°(c2-r2-)V2. 

According to ( p.l3| ) and ( |3.1| ), we need to determine the behaviour of r as proper time r -^ 0. 
(As usual, we fix the origin of proper time so that r(0) = 0.) The previous equation may 
be solved asymptotically by expanding the right hand side and then inverting the resulting 
integral with the result 

r = c|r|^/('"+^)+0(|r|^) (4.2) 
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where j3 > l/(2a + l). 

Then the tangential Jacobi fields have the asymptotic behaviour 

rx ~ c|r|2"/(2"+i) (a ^ 1/2), 

rx ~ c|r|^" In |r| (a = 1/2). 

The behaviour of the radial Jacobi fields is dictated by ( [4. 1| ) which from the above has 
the behaviour 



F(t) ~ cilrl 



-2 



where ci is a negative constant. Then in the notation used above, Go is a negative constant, 
and so this is a type (RSPlc) singularity. The asymptotic behaviour of a is 

where Wi,2 = (1 ± (1 - 4^0)^/^)72, and so 

\\V{t)\\ ~ l/g|r|-2+4a/(2a+l)_ 

Therefore, for any value of a, there will exist radial timelike geodesies along which ||V^(t)|| 
diverges, has zero limit and has finite limit as the singularity is approached. These different 
possibilities arise from the different choices available for c in ( [1.2| ) which give the value 
Go = — 4a(2a + l)c~^"~'^. Starting from some fixed value r = r.,, at r = r.,, < 0, we see that 
c is essentially a measure of the initial velocity of an observer falling radially inwards from 
r*. By tuning this velocity, an observer could in principle ensure that his ||V"(r)|| is finite 
in the approach to the singularity. However, in practice this would be of little help to him 
since as pointed out above, the observer experiences infinite tangential crushing and radial 
stretching in the infall. Furthermore, his initial velocity would have to be tuned with infinite 
precision to obtain 7^ limT-^o- ||^(t)|| < 00. According to the definition above, this is a 
strong singularity. 

D. Marginally bound spherical dust 

The marginally bound LTB space-time (spherically symmetric inhomogeneous dust) has 
line element 

where the prime indicates differentiation with respect to the coordinate 1]. For the collapsing 
case, 

r%r],t) = -m{rj){to{v) - t)^ 



where m, to are arbitrary functions of 77. See |16| for details. The energy density p of the 



dust, which is proportional to the Ricci scalar, is given by 
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4:7ip 



m' 



Thus as well as the central singularity at r = (occurring when t = to (77), there are so- 
called shell-crossing singularities occurring when r' = [|l^. These generally occur before 
the central singularity, at non-zero radius and so are non-central. It has long been believed 
that these scalar curvature singularities are weak. However, it seems that this has only been 



properly established for null geodesies approaching the singularity |§,|n|l- As we have seen 
above, this weakness is completely independent of the structure and nature of the singularity 
apart from the fact that it is non-central. We fill this gap by proving that all radial timelike 
geodesies terminate in a weak singularity. 

According to Proposition Two, the strength of a shell- crossing singularity is governed by 
i ^A^ ). Using (U), we find that 



r,r „ m' m 
2p^f f = 2 

^f^ JUV 9 , ^ Q • 

ty ii lyi lyKj 



The terms m! jr^ and mjr'^ will both be finite in general in the approach to the singularity, 
the former being positive, assuming positive energy density. Thus the behaviour is governed 
by F = 1/r'. We will show that 

lim — = (4.3) 

T^o- r' 

along any radial timelike geodesic approaching the singularity. This shows that the singu- 
larity is type (RSP3a), and is therefore weak by Proposition Two. 
The radial timelike geodesic equations are 

-P + (r') V = -1 (4.4a) 

r'i] + 2r^iT] + r"?)^ = (4.4b) 

i + r'r'^rf = (4.4c) 

where the overdot indicates differentiation with respect to proper time along the geodesic 
and the subscript is differentiation with respect to the global time coordinate t. Along 
each geodesic approaching the singularity, we choose the origin of proper time so that the 
singularity is at r = 0. 



We find that 



r' = LH+ 2t^ 



3 \m to — t^ 
so that at a shell-crossing singularity, 

m' _ 2t'o 



m to — t 

The following terms will enter into our analysis; 
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9 rf' 
< = ^77%7 (4-5) 



3 (to - 1) 



j\2 y, I ^1 1^ 



2 



^"-^-3l^-i^j-i^-(i;^l- («) 



Then 



where here and subsequently, evaluation at zero means in the limit r — > 0^ along a geodesic. 

Generically, ^"[(0) and r"(0) will be non-zero. If this were not the case, there would be 
extra conditions imposed on m and to for all values of r^, which would result in a loss of 
generality. For example, if r^(0) = 0, then tQ(?7) = for all rj. In this case, the space-time 
is homogeneous and isotropic. The condition r"(0) = imposes less severe but nonetheless 
significant restrictions. So we assume henceforth that t'j(O) and r"(0) are non-zero. 

We also need to track the evolution of r and r' along the geodesies. We have 




rtt + r'ri=\ 1 + r't) (4.8) 



We now prove (|4.3|) , which demonstrates the weakness of the singularity. 

Case One:limT-^o- 1^1 < oo|. 

By (|4.4a| ), t(0) = 1. The sign comes from the assumption that the geodesic is future 
directed and the fact that t is a global time coordinate. The past directed case proceeds 
in an identical manner. By ( [4.9|) , (r') will be finite in the limit r -^ 0^. If this limit is 
non-zero, we can apply I'Hopital's rule to r'^/r' to prove ( ^.31 ). The other possibility is that 
(r')(0) = So now assume this to be the case. 

Suppose further that |//(0)| < oo. Then (r'f/)(0) = 0, and so taking the limit of ([4.4bD , 
we have 



= lim {2r[i + r"r]) 
= lim {r[i + (r')), 

T— ►0~ 

which gives ^^(0) = 0, in contradiction of one of our assumptions. So if (r')(0) = 0, then we 
must have |r)(0)| = oo. 

Using I'Hopital's rule twice, we have in this case 



lim — = 


= lim 


r^o- r' 


r^O- 



We calculate 

[r') = r^^t + r^t + ir^ tr] + r rj + r rj . 
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From ([4.4c| ), t(0) = and the terms r^^, r^' and r'" will be finite in the appropriate limit. 
Thus the dominant term is r"fi, giving limT-^o~ V\ = C)0, proving ( [4 .31 ). 

Case Two: liniT-^^o- 1^1 = oo- 

Suppose that |r'?7|(0) < oo. Then by ( [4.4aD , t(0) is finite and so ( |4.9| ) gives |(^')l(0) = °o- 
We then use I'Hopital's rule to prove (^l3| . 

Finally, suppose that |r'?7|(0) = oo. Then by ( [4.4a| ), 



lim t 




Then by 



lim 



\r" ± r'r[\{0) lim |?7| 



This can be finite only if (r"±rV^)(0) = 0. But this limit is generically equal to r"(0), which 
is non-zero, and so we have |(r')|(0) = oo. Again, I'Hopital's rule is used to prove (^^). 

This completes the proof of ( [4.3|) for all radial timelike geodesies and thus demonstrates 
the weakness of the singularity. 



V. CONCLUSIONS 



The central results here are contained in equations ( 2.12 ), (|2.13 ) and ( |2.15|) . These 



provide a set of covariant equations, the asymptotic solutions to which (which require in- 
formation about causal geodesies) determine the strengths of singularities in spherically 
symmetric space-times. The notion of 'strength' is in a slightly modified form to Tipler's 
original definition [^ ; the modification is clearly motivated and is illustrated by the examples 
in Section 3. 

Proposition One demonstrates the important point that the behaviour of null geodesies 
tells us nothing about the strength of a non-central singularity. Also, a null geodesic ap- 
proaching a central singularity terminates in a strong singularity unless r has a finite, non- 
zero limit at the singularity. Proposition Two lists the possible ways in which strong or 
weak singularities may occur. 

In addition to studying the toy models, we were able to demonstrate conclusively the 
weakness of the naked singularity in Roberts' space-time and the shell-crossing singularities 
in collapsing spherical dust. This latter proof shows that while detailed qualitative infor- 
mation about causal geodesies is required, we do not need the full solution of the geodesic 
equations. Therefore there is good hope that the results above may be successfully applied 
to other situations. 

One of these is the case where the singularity occurs at a point where the metric is 
continuous and non-degenerate {det{gab) 7^ 0) (we will refer to such as a continuous non- 
degenerate singularity). It seems plausible that in such a situation, the singularity must 
necessarily be weak. The argument goes roughly as follows. Solutions of the time-like 
geodesies of section 2 typically behave as u = 0(1), v = 0(1) as r —>■ 0. Then to obtain 
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a strong curvature singularity, the Riemann tensor components must diverge faster than 
OIt'"^); integrating twice cannot yield a finite metric. However this argument might not 
hold for an (ISPS) singularity, and perhaps not for other cases. A careful analysis of ( p.6| ) 
and ( p.l3|) should be able to yield either a theorem stating that a continuous non-degenerate 
singularity is indeed weak, or produce examples to the contrary. The statement that a 
continuous non-degenerate singularity is necessarily weak has been made, or the conclusion 
been used, on several occasions in the literature in connection with studies of the Cauchy 
horizon singularity in black holes and singularities in plane wave space-times. This has 
usually been accompanied by separate calculations verifying that the singularity is indeed 



weak |T9|J20|] , but this has not always been the case [pT|-|23||. Thus it appears to be of 
importance to determine exactly when one can conclude weakness for a continuous non- 
degenerate singularity. 

Clarke and Krolak |T^ have given necessary and sufficient conditions, in arbitrary space- 



times, for a singularity to be strong, the conditions involving integrals of certain curvature 
terms along geodesies. An advantage of our work is that it deals with the full set of Jacobi 
fields Jr rather than the volume element V{t). As the toy model of Section 3.3 shows, this 
can be important. Also, the decisive term here 26^-^/™ is slightly simpler than the decisive 
terms in [TO] . It may be possible to use the results here to investigate the connection between 



Tipler's definition of strengths of singularities and Krolak's limiting focusing conditions [22 
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